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By Larry D. Pinson 
Langley Research Center 


SUMMARY 


An analysis using linearized membrane theory has been made to obtain spring 
constants for ellipsoidal bulkheads to be used in longitudinal vibration analy- 
ses of liquid- propellant launch vehicles. A closed-form solution is presented 
for the volume increment and first moment of the volume increment for two types 
of loading, namely, constant pressure and hydrostatic pressure. Plots are pre- 
sented from which the volume increment, first moment of the volume increment, 
and spring constants for ellipsoidal bulkheads can be obtained. Use of the 
equations and plots is illustrated by an example. 


INTRODUCTION 


Stability problems involving interaction between the propulsion system and 
the structure have recently been encountered in launch vehicles. These problems 
have led to increased interest in the accurate determination of longitudinal 
modes and frequencies of the vehicles. However, knowledge of the elastic char- 
acteristics of structural components is required before the modes and frequen- 
cies can be calculated. 

Among the structural components whose elastic properties must be deter- 
mined are the propellant tanks . These tanks, for liquid-propellant launch 
vehicles, are typically cylindrical with ellipsoidal ends. Wood, in reference 1, 
presents a spring-mass representation of such a tank and shows how the spring 
constant can be determined by considering the bulging of the cylindrical part 
of the tank and neglecting the flexibility of the bulkhead. Wood then shows how 
the spring constant for the flexible bulkhead, once obtained, may be incorpo- 
rated into this spring-mass model. The model of reference 1 is reproduced in 
figure 1. The spring constant associated with the bulkhead is determined from 
a consideration of the structural properties of shells of revolution. Treating 
the shell as a membrane, this consideration involves the determination of the 
volume change and first moment of the volume change due to pressure loading. 
Sylvester, in reference 2, derives a general nonlinear expression for the volume 
change in a shell of revolution and applies the linearized form of this expres- 
sion to an ellipsoidal bulkhead subjected to constant internal pressure and 



(a) Geometry and (b) Pressure distribution. (c) Equivalent spring- 

nomenclature. mass model. 

Figure 1.- Ellipsoidal bulkhead under hydrostatic pressure and its equivalent spring- 

mass model. 

hydrostatic internal pressure. Sylvester did not, however, derive expressions 
for the first moment of the volume increment. 

The purpose of the present investigation is to determine the spring con- 
stant for an ellipsoidal bulkhead for use in a spring-mass model such as that 
presented in reference 1. The spring constant for an ellipsoidal bulkhead is 
determined by considering the change in volume, or volume increment, and the 
first moment of the volume increment by using an approach somewhat different 
to that of reference 2. Constant and hydrostatic pressures are considered in 
determining the deformations of the bulkhead, and the results are presented as 
plots which may be used to determine the spring constants for ellipsoidal bulk- 
heads with arbitrary depth- to- radius ratios. 


SYMBOLS 


A cross-sectional area of cylindrical portion of tank 

a length of semimajor axis 

b length of semiminor axis 

C constant of integration 

c positive constant denoting number of g units of axial acceleration 
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E Young’s modulus 

F(n, v) function defined by equation (26) 

G(n, v) function defined by equation (29) 

g acceleration of gravity 

H(n, v) function defined by equation (37) 

h height of liquid in cylindrical portion of tank 

k spring constant 

M first moment of liquid volume 

first moment of volume increment 
mp total mass of propellant 


N 0 

N0 

n 

P 

q. 


stress resultant in circumferential direction 
stress resultant in meridional direction 

depth- to- radius ratio for ellipsoid of revolution, b/a 
arbitrary pressure 

ratio of liquid height h to tank radius a 
r 2 sin 0 


r l 


radius of curvature in meridional direction 


r 2 radius of curvature in circumferential, direction 

S surface 

t bulkhead thickness 

V total liquid volume 

AV volume increment, change in volume of bulkhead due to its deforma- 
tion under load 

v displacement in meridional direction 

W total weight of liquid 

w displacement normal to shell, positive outward 
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z distance measured from base of shell 

unit strain in circumferential direction 

unit strain in meridional direction 

tf distance from free surface of liquid volume to centroid of volume 

before deformation 

change in location of centroid of liquid volume due to deformation 
of bulkhead 

9 coordinate of longitude 

v Poisson's ratio 

| dummy variable 

p mass density 

0 colatitude 

a) natural circular frequency 

Subscripts: 

c constant pressure 

H hydrostatic pressure 

z reference to coordinate z 

£ reference to distance from free surface of liquid volume before 

deformation 


ANALYSIS 


The spring constant for an ellipsoidal bulkhead which is part of a liquid- 
propellant tank is determined by using linearized membrane theory for shells. 
The bulkhead is assumed to be attached to a cylindrical shell as shown in fig- 
ure 1. The depth of the liquid in the cylindrical portion of the tank is h, 
the radius of the cylinder and the length of the semimajor axis of the ellip- 
soid is a, and the length of the semiminor axis of the ellipsoid is b. Only 
the displacements of the ellipsoid are considered herein - the contribution of 
the cylinder to the longitudinal spring constant is treated elsewhere (ref. 1, 
for example). The assumptions used herein are the same as those of reference 1; 
however, a brief discussion of the model is given below. 
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Spring -Mass Model 


The spring-mass model used to represent the propellant -tank combination is 
shown in figure l(c). The deflection of the mass mp is taken to be the change 
in location of the centroid of the liquid due to the tank deformation. The 
liquid is assumed to be incompressible and to act essentially as a rigid body 
insofar as it becomes allowable to assume that all its mass is concentrated at 
its centroid. The cylindrical shell is assumed to be very thin so that bending 
stresses are negligible in calculating deflections. The spring constant k 
represents the stiffness of the bulkhead. When this spring constant is found, 
it is incorporated into the model as a spring in series with another spring 
calculated on the basis of an infinitely rigid tank bottom. The result of this 
combination is the spring kp. On a qualitative basis, the model corresponds 

approximately to the propellant-tank interaction. 

It is the purpose of this paper to establish values of k when the bulk- 
head is an ellipsoid of revolution. 


Spring Constant 

For the purpose of finding the spring constant for the bulkhead, the 
cylindrical portion of the tank is assumed to be infinitely rigid. If the 
liquid is then assumed incompressible, any movement of the centroid will be the 
result of deformation of the bulkhead. 

The loading on the bulkhead consists of a hydrostatic component of pres- 
sure and a constant component equal to the pressure due to the weight of the 
liquid above the bulkhead- cylinder connection. The loading due to tank pres- 
surization (ullage pressure) is assumed to be constant with time and therefore 
may be neglected in the dynamic analysis. 

The spring constant k used herein is defined as the ratio of the total 
force acting on the bulkhead due to the weight of the liquid to the displace- 
ment of the centroid of the liquid volume. This total force is equal to pcgV 
where V is the volume of the mass of liquid, c is a positive constant 
denoting the number of g units of acceleration, and p is the mass density. 

With this definition of k, the problem is to find the resultant force 
acting on the bulkhead and the displacement of the centroid of the mass of the 
liquid due to the deformation of the bulkhead. 

The resultant force is determined as the weight of the mass of liquid in 
the tank under an acceleration relative to the tank. That is 

W = pcgV (1) 

where W is the weight of the mass of liquid. 

The change in centroidal distance is determined by considering the volume 
change of the bulkhead and the first moment of this volume change. 
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The equation which locates the centroid of the liquid is 


where £ is the distance from the free surface of the mass of liquid to its 
centroid before tank deformation and is the first moment of the volume 

occupied by the mass of liquid with the reference axis located at the free sur- 
face of the mass of liquid. (See fig. 1.) 

After the bulkhead deforms, the volume of the bulkhead will change by some 
amount AV. The free surface of the liquid in the tank will drop by an amount 
which gives a volume change AV in the cylindrical portion of the tank. Let 

be the first moment of the volume change in the bulkhead about an axis 
located at the junction of the bulkhead and cylinder. Then the change in the 
first moment of the volume about an axis located at the free surface of the mass 
of liquid before tank deformation is 


AM£ = m z + bAV 


(AV) 2 

2jta 2 


( 3 ) 


where is the change in first moment of the volume occupied by the mass of 

liquid about an axis located at the free surface of the mass of liquid before 
tank deformation, and a is the radius of the cylinder. Since the volume of 
the liquid does not change, the change in centroidal distance is 



(AV) 2 

2jta 2 V 


(*) 


where A£ is the change in the centroidal distance of the mass of liquid due 
to deformation of the bulkhead. 

By the previous definition the spring constant for the bulkhead is 


k 


pcgV 2 


m z + hAV - 


(AV) 2 

2jra 2 


( 5 ) 


where k is the spring constant. The quantity (AV) / 2jta 2 is small when com- 
pared with the other terms of the denominator of equation ( 5 ) and may be 
neglected so that 

_ pcgV 2 
AMg + hAV 


( 6 ) 
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The neglect of this quantity is consistent with other linearizing assumptions 
made herein. 

The quantity V is found from the geometry of the tank but the quantities 
AM Z and AV must be found by integrating the normal displacements due to the 
load over the surface of the bulkhead. 

Volume Increment 

The equation for the volume increment AV for a membrane shell of revolu- 
tion with arbitrarily shaped meridian is derived in this section. A typical 
shell of revolution and the coordinate system is illustrated in figure 2. 

It is assumed in the following derivation that bending stresses are 
negligible and that changes in the radii of curvature are negligible. A con- 
sequence of the assumption that bending stresses are negligible is that there 
is no control over the normal displacement at the boundary of the shell. This 
condition makes it impossible to match the displacements of the bulkhead and 
the cylinder at the boundary. It is considered however that this is a localized 
effect and will not substantially affect the overall distribution of the 
displacements . 



Figure 2.- Typical shell of revolution- 
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From figure 2 the volume increment due to a loading which produces a nor- 
mal displacement w is 


AV = 2ix / r o r}Wd0 

J h 


where 


r 0 - r 2 sln t> 


r^ is the principal radius of curvature in the meridional direction 0; T2 

is the principal radius of curvature in the circumferential direction 0j and 
02 and 02 are the lower and upper limits of the shell, respectively. Refer- 
ence 3 gives the displacements for axi symmetric loading as: 

w = r 2 € Q - v cot 0 (9) 


..in 4 f ^ 

r \J sii 


d0 + C 


where C is a constant of integration which is determined from the requirement 
that the meridional displacement v have some known value at a particular 
value of 0 and 


*( 0 ) - 


i\ew - r 0 e. 


In these equations ana €q denote unit strains m tne. jienaionai ana cir- 

cumferential directions, respectively. Substitution of the equations for r Q , 

w, and v into equation (7) yields the following expression for the volume 
increment 


denote unit strains in the. meridional and cir- 


rh 2 w w 

= 2n / r l r 2 € 0 sin 0 10 “ 2n 

J 01 


sin 0 


d0 + C sin 0 cos 0 d0 


If it is required that the meridional displacement v be zero at the lower 
limit of the shell, the constant C may be eliminated by putting limits of 
integration on the bracketed integral. Introduce a dummy variable of integra- 
tion into the integrand of the bracketed integral and place the limits 0p and 

0 on the integral. The result is 
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sin 0 cos 0 d0 


AV 


St 



v •„ 2 

r l r 2 e 0 


sin 0 d0 - 2jt 



f(S) 

sin | 



where £ 


is a dummy variable of integration. 


( 13 ) 


First Moment of the Volume Increment 

In the development of the equation for the first moment of the volume 
increment, the same basic assumptions -which were made for the development of 
the equation for the volume increment are made. 

From figure 2, the first moment of the volume increment about the base of 
the shell due to a loading which produces a deflection w is 

o02 

AM = 2jt / r o r l zw d0 (l^) 

J fa 


where z is measured from the base of the shell. Following the same procedure 
as used in the previous section for determining volume increment, the expres- 
sion for the first moment of the volume increment about the base of the shell 
for axisymmetric loading is found to be 


AM = 2n 


pfa ? . . Pf 2 

r l r 2 ze 9 sln P d 0 " 2jt L r l 

J fa J fa 


,02 


r^z 


ft f(6) 

0 X sln 6 


dl 


sin 0 


cos 


0 d 0 

(15) 


where the requirement has again been made that the meridional displacement v 
is zero at the base of the shell. 


Application to the Ellipsoidal Bulkhead 

In order to apply the general expressions for the volume increment and the 
first moment of the volume increment (eqs. ( 13 ) and (l 5 ))> several quantities 
which are interrelated through the loading and shell geometry must be estab- 
lished. The quantities in equations ( 13 ) and ( 15 ) which are related to the 
loading are the unit strains and . These strains are related to the 

stress resultants through Hooke’s law and the stress resultants are related to 
the loading through the equilibrium of forces on the shell and the geometry of 
the shell. Specifically, Hooke’s law states that 


( 16 ) 


e i> - St (■¥ - vN e) 
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and 


€ ® = Et ( Ne " VN ^) <17) 

where E is Young's modulus; t is the thickness of the shell; v is Poisson's 
ratio; and and Nq are the stress resultants in the meridional and cir- 

cumferential directions, respectively. To find the stress resultants the radii 
of curvature of the shell must be known. Also needed to apply the expression 
for the first moment of the volume increment is the quantity z which is the 
distance from the base of the shell to any point on the bulkhead. For the 
ellipsoid of revolution these quantities are: (See ref. 3 * ) 


n^a 

(sin 2 0 + n 2 cos 2 0 

) 3/2 

a 


^sin 2 0 + n 2 cos 2 0 

r 

-n 2 a cos 0 


fsin 2 0 + n 2 cos 2 0 

\l /2 


where n is the depth- to- radius ratio of the ellipsoidal bulkhead. The 
pertinent range of 0 is ^ % 0 it . 


(18) 


( 19 ) 


( 20 ) 


Stresses in an ellipsoidal bulkhead .- The stress resultants for the 
ellipsoidal bulkhead are found from the equilibrium equations and the external 
forces. These forces are related to the pressure distribution. The loadings, 
or pressures, being considered herein are divided into two components as illus- 
trated in figure 1 - a constant normal pressure equal to pcgh and a hydro- 
static pressure equal to pcgz. For these loadings, the membrane theory for 
shells of revolution (ref. 3 , for instance) gives the following stress 
resultants for the two components of the load on the ellipsoidal bulkhead: 

For constant pressures - 


N0 = 


pcgha 

2 ^sin 2 0 + n^cos^0^ ^ 


( 21 ) 


and 


N 0 



pcgha 


p P 
+ n^cos c ' 




1/2 


pcgha /~ s j _ n 20 + n 2 cos 2 ^j 
2 n 2 


( 22 ) 
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For hydrostatic pressure - 


and 


pcga 2 n 

n^cos-^0 + 1 

(si + n^cos^0 

f 

3 

sin20 (sin^0 + n 2 cos^0) 



(23) 


-pcga 2 n 2 cos 0 

pcga 2 

3 /2 

n^cos-^0 + (sin^0 + n^cos 2 0) 

( sin^0 + n 2 cos 2 0) 

3 n 

_ sin^0 


( 24 ) 


Having obtained the stress resultants, the strains may be found by sub- 
stituting the stress resultants for the particular component of loading being 
considered into equations (l 6) and (17) • 

Volume increment for an ellipsoidal bulkhead .- The volume increment for an 
ellipsoidal bulkhead is obtained by substituting the strains for the appropriate 
loading condition into equation ( 13 ), the general expression for the volume 
increment. If the thickness of the bulkhead is variable the strains must con- 
tain t as a function of 0 . In the present analysis the bulkhead thickness 
is assumed to be constant. It is also assumed that the bulkhead is one- half 

of an ellipsoid of revolution, thus giving the range of integration 75 = 0 = rt . 

The integrations are straightforward, and give the following results for the 
volume increment. 


For constant internal pressure, after making the substitution £ = cos 0 , 
it is found that 


AV C = 


itpcgha 

2 Et 


|i(n,vj] 


( 25 ) 


where 


In 


r l + \]l 


n 


F(n,v) = 


15 - 12v 

~E 


_ 3 _ 

4 n £ 


V 

(15 " 12v)n 2 0 + 1 

1 

V 

r 

L 8 V n 2 J 

1 2 { 


- yi - n c 


;^1 - n £ 


(26) 


The function F(n,v) is plotted against n in figure 3 for various 
values of Poisson’s ratio. Equation (25) agrees with the result of reference 2 . 
The volume increment for a hemispherical dome can be found by taking the limit 
of AV C given in equation (25) as n approaches unity. The result for a 
hemisphere is found to be 


AVc = - V) 


(27) 


11 



The simplicity of the geometry of the hemisphere makes the calculation of 
the volume increment for this case quite simple and the calculation provides an 
independent check, though not an absolute proof, of the validity of the expres- 
sion for the volume increment for the ellipsoidal bulkhead. An entirely sepa- 
rate calculation was made for the hemisphere and the result agreed with 
equation (27). 
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Figure 3.- Variation of F(n,v) with depth- to- radius ratio for various values of 

Poisson’s ratio. 









For hydrostatic pressure, after making the substitution | = sin it is 
found that 


AV H = 


npcga- 

2Et 


[d(n,v)] 


where 


( 28 ) 


G(n,v) = 


-13 + 37 n 2 + 10n 5 + f~— 

1 + n 2vn j 

9 


4 + 


n^- 
1 + n. 


+ 



(29) 


Equation ( 29 ) is plotted in figure 4 against n for several values of 
Poisson^ ratio. The volume increment for the hemisphere subjected to hydro- 
static pressure is found by letting n approach unity in equation ( 28 ). The 
result is 


- V) (50) 

Separate calculation of the volume increment for the hemisphere subjected to 
hydrostatic pressure yields a result which is in agreement with equation ( 30 ). 

Equation ( 29 ) differs from the result presented in reference 2 in that 

P P 

the quantity 2vn in the last term of the equation is vn in reference 2. 

Communication with the author of reference 2 confirmed the validity of the 

expression as given in equation ( 29 ). 

First moment of volume increment fo r an ellipsoidal bulkhead .- Before 
proceeding with the development of the relations for the first moment of the 
volume increment the following theorem is considered which is presented in ref- 
erence 4 on the first moment of the volume increment in any shell: Let p be 

an arbitrary pressure distribution and wg be the normal displacements due to 
a hydrostatic pressure resulting from a liquid whose unit weight is peg. The 
first moment of the volume increment due to the arbitrary pressure distribution 
is the surface integral 


AM = 



(31) 


This theorem can be used to determine the first moment due to constant 
pressure. If p is selected to be the constant pressure pegh, equation ( 31 ) 
becomes 

= h I v H dS (32) 

J S 
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Therefore, the first moment of the volume increment due to constant pressure is 

m c = hAv H (3*0 

Substitution of equation (28) into equation (34) gives the first moment of the 
volume increment created by the constant component of pressure 


«<c - 2fi^[a(n,v)] (J5) 

where G(n, v) is given by equation ( 29 ). 

Unfortunately no such simple relationship exists for the first moment of 
the volume increment caused by the hydrostatic component of the pressure and 
one must resort to the method used in obtaining the volume increment. 


To find the first moment of the volume increment which is caused by the 
hydrostatic component of the pressure, the equations for the stress resultants, 
equations ( 23 ) and (24) are substituted into equations (l6) and (17) to obtain 
the strains, and the resulting expressions are substituted into equation ( 15 ). 

This expression is then integrated over the range ^ s 0 tt and it is found 
that 


m E = 


rrtpcga^ 

2Et 


|H(n,vj] 


(36) 


where 

H(n,v) 


^[27 - l 6 n - v(2kn - 3 )] 


2 - 6n 2 + 3n 4 ' - vn 2 (2 - n 2 ) 


In 


ll + \Jl - n 2 


\jl - n 2 


9^1 - n 2 


n (3 - v) 
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2]jl - n 2 


n 2 lnl 


+ Vi 


n c 


- \Jl - n 2 


2(l - n 2 ) }jl _ n 2 


8n (l + v) /1 + n 
+ r In 

n 


( 37 ) 


The function H(n, v) is plotted in figure 5 against n for various 
values of Poisson^ ratio. 

When equation ( 36 ) is specialized to the hemisphere by letting n approach 
unity, the following expression results for the first moment of the volume 
increment in a hemisphere subjected to hydrostatic pressure: 
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Figure 5.- Variation of H(n,v) with depth- to- radius ratio for various values of 

Poisson’s ratio. 






•( 38 ) 


^ + sup) ln(2 ^| 

Separate calculations for the hemisphere, as in the checks on the volume 
increment, yield an expression which is in agreement with equation (38). 


Spring Constants for Ellipsoidal Bulkheads 

De velopment of equations .- With expressions for the volume increment and 
first moment of the volume increment, the equation for the spring constant can 
be developed according to equation (6). 

The volume of liquid within the tank is 

V = ^ 3<1 + 2n) (39) 


where q is the ratio of the height of liquid within the tank h to the radius 
of the tank a, and n is the depth- to- radius ratio of the ellipsoidal 
bulkhead. 


Adding the volume changes for the two components of loading and the first 
moments of the volume changes and substituting these sums, along with equa- 
tion (39) f or the volume into equation (6) yields the following expression for 
the spring constant for the ellipsoidal bulkhead 


k = Et - 


<3U;>q 


9 (h(u, v) + 2 qG(n,v) + q 2 F(n,vi 


( 40 ) 


where the functions F(n, v), G(n, v), and H(n, v) are defined by equations ( 26 ), 
(29), and (37)> respectively. 

Equation ( 40 ) gives the spring constant for an ellipsoidal bulkhead sub- 
jected to constant and hydrostatic internal pressure. This equation shows that 
the spring constant is a function of the geometry of the ellipsoid - the 
ratio nj the height of the liquid in the tank - the ratio qj Young's modulus ; 
the bulkhead thicknessj and Poisson's ratio. The spring constant does not 
depend on the radius of the tank directly - the radius affects the spring con- 
stant only through the ratios n and q. 


In certain launch vehicles the oxidizer and fuel are separated by a bulk- 
head which is common to both tanks, and geometry similar to that shown in fig- 
ure 6 results . The possibility of buckling due to the resulting compressive 
stresses in the bulkhead is offset by pressurizing the aft tank. For this con- 
figuration the volume change and first moment of the volume change may be cal- 
culated by using the functions already developed by finding the volume change 
or first moment of the volume change due to a constant pressure pcgh and 
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subtracting the volume 
change or first moment of 
the volume change due to 
a hydrostatic pressure 
(shown by the dotted line 
on the pressure distribu- 
tion of fig. 6). 

The equation for the 
volume of the liquid in 
the tank is 


V = 2|^(3q - 2n) (1+1) 

where the quantities V, 
a, q, and n have the 
same meaning as before. 

Repeating the pro- 
cedure used in the devel- 
opment of equation (40) 
and noting that the terms 
associated with the hydro- 
static component of pres- 
sure are subtracted, 
rather than summed, the 
following equation for the 
spring constant results: 




Pcgb 


(a) Geometry for 
an inverted 
bulkhead. 


(b) Pressure 
distribu- 
tion. 


Figure 6.- Common ellipsoidal bulkhead under hydrostatic 
pressure. 


k 


Et 


2n(3g - 2nf 

9[H(n,v) - 2qG(n,v) + q 2 F(n,v)] 


(42) 


The range of q in this case must be restricted to the range q ^ n because 
it is assumed that the liquid completely covers the bulkhead. 

If the height of liquid is very large compared with the depth of the 
ellipsoidal bulkhead, the contribution of the hydrostatic component of pres- 
sure to the spring constant will become small compared with the constant com- 
ponent of pressure and the expression for the spring constant will approach an 
aysmptotic value. This asymptotic value is found by taking the limit of the 
expression for the spring constant as q increases without bound. Taking this 
limit yields 


k 


= Et 


at 

F(n,v) 


( 43 ) 


Properties of equations .- If equation (40) is divided on both sides by the 
quantity Et, differentiated with respect to q, and the derivative set equal 
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to zero, an equation results which gives the value of q for which the quantity 
k/Et is a maximum. This equation is 


3 H(n,v) - 2 nG(n, v) 
3 G(n,v) - 2 nF(n, v) 


(W) 


If this value of q is plotted against n for a value of Poisson 's ratio of 
l/3, two branches result as shown in figure 7 - These branches become asymptotic 
to the value n = 0 . 827. A similar procedure applied to equation ( 42 ) yields 
the same expression except the sign of the right side of the equation is changed 
from negative to positive. The branch having negative values is of no interest 
in the discussion of equation ( 40 ) because the analysis does not hold when the 
liquid surface is below the ellipsoid-cylinder connection in which case a sepa- 
rate problem must be solved. How- 




ever, for values of n greater 
than 0.827, positive values of q 
exist where the quantity k/Et 
is maximum. Figure 8 shows this 
effect. The curve for n = 1 , 
that is, the hemisphere has a 
maximum at q = 1 . 48 . For values 
of n less than 0.827, no rela- 
tive maximum exists. Thus, for 
values of n less than 0.827, 
and for a Poisson 1 s ratio of 1/3, 
the spring constant will not 
exceed the asymptotic value given 
by equation ( 43 ). In discussing 
equation ( 42 ), opposite conclu- 
sions may be drawn. Since the 
sign of the right side of equa- 
tion ( 44 ) is changed the negative 
branch of the plot of figure 7 
becomes the positive branch and 
the positive branch becomes nega- 
tive. Thus, for values of n 
greater than 0.827, no finite 
value of q exists for which the 
quantity k/Et is a maximum and 
it will never exceed the asymp- 
totic value given by equa- 
tion ( 43 ) . For values of n 
less than O.827, positive values 
of q which are greater than n 


J exist such that the quantity 
k/Et is maximum. This effect 
is shown in figure 9* 


Figure 7*- Variation of critical value of height- 
to- radius ratio vith depth- to- radius ratio for 
a Poisson's ratio of l/3* 


The derivation did not 
include the deformations of the 
cylindrical portion of the tank. 
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Figure 8.- Variation of spring constant for ellipsoidal bulkhead with height-to- radius 
ratio for several values of n and Poisson's ratio of l/3* 

Illustrative Example 


In order to clarify the procedure used in obtaining the spring constant 
and the process of combining the bulkhead spring constant with that of the 
cylindrical portion of the tank, an exanqple is presented using typical 
parameters . 

Suppose it is desired to find the equivalent spring-mass system for a 
tank having the following parameters: 
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The quantity q is ^ = O. 75 , From figures 3 , 4, and 5, it is found that 

60 

F(n,v) = 1.150, G(n,v) = 0.695 j and H(n, v) = O. 615 . Substitution of these 
values into equation (40) yields a spring constant for the bulkhead of 
2.034 X 10 6 lb/in. 


The spring constant just determined is only part of the complete stiffness 
representation for the tank-liquid combination. The complete representation, 
from reference 1, is illustrated in figure 1. The equations used to obtain the 
various spring constants are taken directly from reference_l; however, equa- 

vKp K2 

tion A. 21 on page 8 l of the reference should read Kb = , not Kq = 

1 - v 1 - v 


k l = 


OJ 
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0J 

1 AE 

3 - 2v 2 y 

1 h 


2 - 2v \AE , ,, 

2 l k 


ko = 


3 - 2 v 

vkq 

1 - v 


k 3 = f ~ vk l 


(45) 

(46) 

(47) 


Substituting the given values into equations (45), (46), and (47) and assuming 
the cylindrical tank wall to have the same thickness as that of the bulkhead, 

the following results are obtained: kq = 1.011 x 10 ^ lb/in., 

k 2 = 0.506 X IQ 6 lb/in., and kj = 3.852 x 10^ lb/in. 


The effect of incorporating the spring constant for the bulkhead into the 
spring-mass model of figure 1 may be evaluated by comparing the values obtained 
from equations (45), (46), and (47) with those obtained by assuming that the 
bulkhead is infinitely stiff. If the bulkhead is rigid equations (46) and (47) 
are valid but equation (45) becomes 


_ (2 - 2v) AE 

1 ‘ (T-iv £ ) ^ 


(48) 


By using equation (48) along with equations (46) and (47), the following results 
are obtained for the system with a rigid bulkhead: kq = 2.011 x 10 ^ lb/in., 
k 2 = 1.005 X 10^ lb/in., and k^ = 3-519 X 10 ^ lb/in. The effect of including 
bulkhead flexibility, therefore, is a reduction of 50 percent in both kq and 
k 2 and an increase of 9 percent in k^. 


The effect of bulkhead flexibility on the frequency may be determined by 
assuming the ends to be fixed. The uncoupled liquid frequency is 
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where mp is the actual liquid mass. Since the frequency varies with the 
l /2 power of the sum kj + k 2 , the effect of including the bulkhead spring con- 
stant is a 42-percent reduction in the frequency for this example. 

CONCLUDING REMARKS 


An analysis was made to obtain spring constants for ellipsoidal bulkheads 
to be used in longitudinal vibration analyses of liquid-propellant launch vehi- 
cles. The analysis was based on linear membrane theory of shells- Closed-form 
expressions for the spring constant for ellipsoidal bulkheads subjected to con- 
stant and hydrostatic pressure have been presented. The volume increment and 
first moment of the volume increment for these pressures were derived. Plots 
are presented to aid in the evaluation of the spring constants for bulkheads of 
arbitrary depth- to- radius ratios. The relation of the spring constant of the 
bulkhead to the entire tank model has been discussed. The procedure has been 
illustrated by means of an example. The example showed that including the bulk- 
head flexibility caused a significant change in the frequency of the particular 
system considered. 

The work presented herein could be extended by computing the spring con- 
stant for other cases, among which are the following: (l) the case in which 

the level of the liquid is below the ellipsoid- cylinder connection, or ( 2 ) by 
computing the spring constant for another configuration, say a conical shell. 

The work could further be extended by directing efforts toward the improvement 
of the propellant-tank model. 


Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 13 , 1964. 
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